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ABSTRACT 



Aims. Hot Jupiters are thought to belong to single-planet systems. Somewhat surprisingly, about one quarter of hot Jupiters were 
reported to exhibit Transit Timing Variations (TTVs). The aim of this paper is to identify the origin of this observation. 
Methods. We present TTV frequencies and amplitudes of hot Jupiters in Kepler QO-6 data with Fourier analysis. 
Results. We identified 29 systems with TTV above 4cr confidence, about half of them exhibiting multiple TTV frequencies. Thirteen 
of these objects (HAT-P-7b, KOI-13, 127, 188, 196, 203, 225, 254, 428, 607, 774, 897, 1176) likely show TTVs due to a systematic 
observational effect: long cadence data sampling is regularly shifted transit-by-transit, interacting with the transit light curves and 
introducing a periodic bias. In case of other systems, the activity and rotation of the host star can modulate light curves and explain 
the observed TTVs. Excluding the systems that were inadequately sampled or showed signs of stellar rotation, we ended up with 
16 systems. Nine of them show only one significant period (KOI-131, 186, 256, 882, 897, 1003, 1152, 1448, 1543), others exhibit 
multiple periods (KOI-823, 977, 1285, 1382, 1452, 1540, 1546). The number of significant periods ranges from 2 (KOI-977, 1382) to 
10(KOI-1546). 

Conclusions. The few hot Jupiters with periodic TTVs that we cannot explain with systematics from observation, stellar rotation, 
activity or inadequate sampling may be good candidates for hot Jupiters in multiple systems or even hosting an exomoon. 

Key words, planetary systems - stars: binaries: eclipsing - techniques: photometric 



1. Introduction 

Transit Timing Variation (TTV) is a major diagnostics of various 
system parameters of extrasolar planets. In multi -planet systems, 
planets perturb each other, leading to correlated TTVs of them 
(Holman et al. 120101 Lissauer et al. 1201 11 1. TTVs have also un- 
covered the presence of further non-transiting planets in solar 
systems (Ballard et al. l20Tn Ford et al. I2012bl l. 

According to our current view, hot Jupiters occur as single 
planets, since they have been not detected in multi-planet sys- 
tems. This picture suggests that hot Jupiters occupy unperturbed 
orbits, hence their orbital motion is Keplerian, and they exhibit 
strictly periodic transit times. Contrary to this picture, current lit- 
erature reports a considerable number of hot Jupiters with TTV, 
which are often periodic (StefFen et aI. I2012a[ Ford et aI. I2012al) . 

In this work we publish TTV periods, TTV amplitudes and 
significance levels for hot Jupiter candidates in Kepler data. For 
the analysis, transit times covering Q0-Q6 were taken from Ford 
et al. |22l2£C] (see Steff^en et al. I2012bl for more details). We 
briefly introduce the most exciting systems and discuss the pos- 
sible source of TTVs for these planets. 



2. Data selection and analysis 

We have selected all single companions from the second list of 
Kepler planet candidates (Batalha et al. I2012I I. and filtered the 
list for planets larger than 6 Re. Then all systems with incom- 



plete data records, too few observed minima (we required more 
than 12 transits of a given system for our analysis) or inaccurate 
orbital periods were neglected. In total, 159 candidates were se- 
lected for further analysis. The median planet size and orbital pe- 
riods were 1 1.23 Re and 5.70 days in the sample, the longest or- 
bital period was 54 days. Thus, the sample contains really close- 
in planets, i.e. they are hot Jupiters. The sample was trimmed fur- 
ther, since some of our targets were reported to be triple systems 
or blends, thus not hot Jupiters (KOI-190 and KOI-609 Santerne 
et al. I2012I I. or the periodic signal was caused by the influence 
of stellar rotation and/or activity (see in Sec. 14.31 ). Three sys- 
tems (KOI- 256, KOI-882, KOI-1540) were labefled in Batalha 
et al. 120121 as showing V-shaped transits. Although planets can 
cause V-shaped transit light curves, but the probability of being 
a false positive is much higher for these objects. Despite this, 
we follow Batalha et al. 120121 in retaining these candidates. Our 
sample is listed in Table [T] The table contains stellar and plane- 
tary parameters taken from the PlanetQuest websit^ These are 
compatible with the values from Borucki et al. 120111 Batalha et 
al. lSbni and the Kepler Input Catalog (Brown et a l.lSOTlT l. Our 
table also shows the vetting flag from Borucki et al. l201 II where 
(2) represents "strong probability candidate, clearly passes tests 
that were applied", (3) stands for "moderate probability candi- 
date, not all tests cleanly passed, but no definite test failure" and 
finally (4) means "insufficient follow-up to perform full suit of 
vetting test." 



http://www.astro.ufl.edu/~eford/data/kepler/ 



http://planetquest.jpl.nasa.gov/kepler/ 
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Fig. 1. Illustration of the sub-Nyquist detection of a super-Nyquist modulation in the orbital motion with Fourier analysis. Top panel; 
one quarter (90 d) long detail from the 600 d long simulations. Dots represent sampling by transits. Bottom panel: Fourier series 
of the 600 d long data set, extended to beyond the Nyquist frequency. The actual frequency and the sub-Nyquist detections are 
highlighted. This specific model had an orbital (sampling) and modulation period of 11. 21 and 4.76 days (which were taken from 
the ZIP code and altitude-above-sea-level of the Budapest station of Konkoly Observatory, as uncorrected random numbers). 



Transit times of these systems were taken from the TTV cat- 
alog of Ford et al. I2012al and were processed by the standard 
Fourier analysis with the MUFRAN package (Kollath 199(3). Here 
we wish to illustrate the most serious drawback of the applica- 
tion of the Fourier approach in such an analysis. Once the orbital 
phase is modulated by a periodic signal, a periodic TTV will 
be observed. That is, a periodic TTV confidently proves a peri- 
odic modulation of the orbital motion. However, the period of the 
modulation cannot be reconstructed from the TTV if it is shorter 
than 2 orbital periods. This is because transits undersample such 
frequent variations, and the observed frequency will be an alias 
of the super-Nyquist peak. This is illustrated in Fig. [T] In the 
general case, this can lead to an order-of-magnitude difference 
between the detected and the actual period. This is unfavourable 
in the sense that one cannot tell practically anything on the pe- 
riod of the TTV, the detected period will signify only a lower 
freqeuncy limit for a periodic process. However, the amplitude 
of the modulation can be properly reconstructed from TTV data 
(except strictly resonant cases), and the TTV amplitude can be 
involved into further analysis. 

The Fourier-transforms were evaluated individually. The 
presence of TTVs were characterized by the detected period, the 
TTV amplitude and the height of the detected peak above the 
noise level. We considered the planet candidates with at least one 
peak higher than 4cr above the grass level as systems exhibiting 
periodic TTVs. 

2A. Virtual frequencies due to regularly-spaced sampling 

We found that in some cases, TTVs exhibited prominent fre- 
quencies that can be easily explained by sampling effects. Most 
of the data are long cadence, and they have a sampling rate of 
29.424 minutes. This is in the order of the duration of one whole 
transit, resulting in 3-6 points belonging to most of the transits. 
Since the orbital period is not an exact multiple of the cadence 
rate, the consecutive transits are sampled at different positions, 
but the lags will evolve very regularly, due to the even sampling 
rate. The instantly realized sampling can introduce biases in the 
transit time determination, and because of the regular evolution 
of the entire process, a virtual variation of transit times may be 




C-frac(P/C) 



Fig. 2. Calculating the virtual frequency in transit timings, due to 
the uneven sampling of consecutive transits. P and C are the or- 
bital period and the cadence, respectively, frac means the frac- 
tional part. 



measured. (A similar feature is seen e.g. in trailed CCD images 
of asteroids, where the very regular - paced - sampling results 
in a virtually undulating motion with a sub-pixel amplitude.) 

In Fig. 2 we illustrate the derivation of this frequency. 
Because the cadence of exposures is (obviously) not synchro- 
nised to transit mid-times, a systematic lag will develop between 
the sampling structure of the consecutive transits. The sampling 
comb will be off one transit later by a time lag of \P - nC\, where 
P and C are the orbital period and the cadence rate, and n is an 
appropriate integer We know that « = [f'/C]^,', where [] means 
the integer part. There are two cases, one when nC is slightly 
less than P (by less then P/2, see this case in Fig. 2), and an- 
other one, when nC is slightly larger than P. In the two cases, 
denoting the fractional part with ] [ brackets, we can write as 
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monoperiodic modulations. Labeled levels show the confidence limits of peak detec- 
dashed lines. In case of KOI-1003 after prewhitening with the 277d peak (f =0.003601 
') does not reach the Act level. 



and following the same logic and inverting the sign because the 
quantity in absolute value brackets is negative 
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Thus, the time lag can be written as C ,s if we define s — 
min(]P/C[, l-]P/C[). The C ■ s time lag is s times one cadence, 
and after l/s occurrences (i.e transits), the initial cadence-to- 
transit configuration will be repeated. This means that the vir- 
tual peak has P/s period, and hence, s/P frequency. Because s 
is always less than 0.5 by definition, the virtual peak will always 
emerge, but will not always be high enough to be detected. In 
summary, if a transiting system exhibits a TTV period near P/ s, 
this period should be considered to be suspicious, because it may 
refer to the virtual period of the system. 

Among other systems, we found that the most prominent ex- 
amples of systems with likely virtual TTV periods are KOI-2 
(HAT-P-7b), KOI-13, and KOI-127. These systems have a TTV 
period peak at 21.3, 5.7, and 25.2 days, respectively, while their 
virtual period is 22.2, 5.7, and 24.3 days. Because we cannot ex- 
clude that the detected TTV period is the virtual one, we omitted 
these objects from our table. 

3. Results 

Twenty-nine systems out of the 159 were identified showing pe- 
riodic TTVs, which is 18% of all candidates. Of these targets, 
17 candidates show TTV peaks convincingly far from the virtual 
periods, and also they were not suggested to be blended objects. 
One more system (KOI-883) was excluded because the TTV pe- 
riod was found to agree with the rotation period of the host star. 
In Figs. [3] &|4l, we plot a sample of Fourier spectra of the found 
TTV variations. 

Surprisingly, 7 of the finally accepted systems exhibited mul- 
tiple frequency peaks (even 9 and 1 3 formally above the 4cr level 
in case of KOI-1452 and KOI-1546, respectively), suggesting 
that the orbital motion is modulated by several periodic pro- 
cesses. The median period of the detected TTV was 50 days. 



that means that at least one half of the planet candidates suffer 
periodic modulations with shorter periods than 50 days. Since 
the detected TTV period is only an upper limit for the actual pe- 
riod of the modulations, very well may be that the median of the 
actual periods is much smaller than 50 days, even an order-of- 
magnitude difference can easily be suspected; but by no ways 
can the median period exceed 50 days. In some cases (e.g. KOI- 
1452 and 1540), the major peak is of very low frequency, sug- 
gesting a long period process. However, these peaks may also 
be aliases of rapid oscillations near some resonance with the or- 
bital motion, thus near to the Nyquist frequency multiplied by 
an even number, and having reflected back near the origin in the 
sub-Nyquist range. The median amplitude of TTVs is 54 sec- 
onds. 

We compared system parameters of hot Jupiters with and 
without TTVs, and found no significant difference in the dis- 
tributions in the 2-dimensional subspaces of the orbital period, 
planet size, stellar mass, stellar radius, semi-major axis. One sig- 
nificant difference was found, but very probably due to numer- 
ical rather than physical reasons: the orbital period of systems 
with periodic TTVs is preferentially shorter than that of all hot 
Jupiters. However, the orbital period is severely biased toward 
the short period range in the selection process, because a se- 
cure detection of any periodic TTV requires many transits, and 
prefers the short orbital periods. 

Since almost all our targets were examined in Steffen et 
al. I2012bl it is worth comparing our results. We note that our 
method differs from that of Steffen et al. I2012bl we fit only a 
single-component sinus at a time, and the definition of the sig- 
nificance level is also different. Instead of taking the uncertainty 
of the amplitude value, which is a local feature, we compare 
the amplitude of the signal to the noise computed for the total 
Fourier-spectrum. We found that the p-values taken from Steffen 
et al. l2012bl are not particularly low for our targets, which means 
that the TTV values of solitary hot Jupiters are not the highest 
in the total Kepler sample. The SNR values agree for about half 
of the common objects, while the significance found with our 
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method is usually higher for the rest of our targets. We attribute 
this difference to the different methods applied. 

3.1. Remarks on individual systems 

Most of our candidates in Table. [T] are listed in Ford et al. 1201 II 
We only refer to that paper in the table if it provides relevant 
information in the context of this work. 

KOI-131 Santerne et al. l2012l observed this object as part of their 
large spectroscopic follow-up program of Kepler planet candi- 
dates. They found a relatively fast rotating exoplanet host star 
with vsini = 27+1 kms"' and place an upper limit of 14.3 Mj 
for the mass of the companion. The authors note that both the 
planetary and blend scenarios are compatible with their data. In 
addition, the object is listed in the Kepler eclipsing binary cata- 
log (Prsa et al. 2012) with twice the orbital period. This is pre- 
sumably because of the slightly different depth of the odd and 
even minima based on the phased light curve on the websit^ In 
(Borucki et al. 1201 11 1 it got a vetting flag 3. KOI-131 is listed as 
a planet candidate in Batalha et al. 120121 as well. We conclude 
that this object is not a strong hot Jupiter candidate. 
KOI-186 This Jupiter-sized planet orbits a solar-type star 
(M=1.06 Mq). We found two significant periodicities, namely 
13.877 and 7.243 days. The later is very close to the presumed 
stellar rotation period (7.84 d), so we are left with a TTV candi- 
date exhibiting a single periodic variation. 
KOI-256 According to Muirhead et al. 120121 the host star is a 
metal-rich M3 dwarf. Their measurement significantly reduced 
the size of the planet candidate from 25.34 to 5.60 Re- The com- 
panion shows a single 41.8 d periodicity in its transit timing. 
KOI-412 Data from Q5 is missing because the target fell on the 
failed Kepler Module 3. The highest peak in the Fourier spec- 
trum of the transit timing data of this object is located exactly at 
the Nyquist-frequency, i.e. twice the orbital period of the planet. 
We can not exclude a signal resulting from a resonant object with 
1 :2 mean motion resonance in the system. The amplitude of the 
TTV is uncertain for the same reason. 



' http://keplerebs.villanova.edu/ 



KOI-822 For similar reasons as in the case of KOI-412, Q5 data 
are is missing. Strong rotational modulation is seen in the light 
curve with about half the detected TTV period, as well as its first 
harmonic, so it may influence the transit timing measurements. 
Therefore we flagged this star as an uncertain case. 
KOI-823 This object fell on the bad module in Q6. It shows 
complicated, multi-peaked Fourier-spectrum. While some of the 
peaks are close to peaks seen in the light curve itself (142.9, 
59.5, 19.0 d), others are found to be independent of stellar and 
sampling effects. 

KOI-882 Single-periodic TTV candidate. The long period seen 
in the TTV (41.9 d) is unlikely to be disturbed by other relatively 
long-period periodicities found in the light curve, the closest one 
is 37.1 d. 

KOI-883 Based on its KIC-parameters, the host star is an early 
type K dwarf. We found two peaks in the Fourier-spectrum of the 
transit times of its companion (6.6 c/"' and 9.1 c/"'). The 6.6 c/"' 
peak is uncomfortably close to the alias caused by the sampling 
effect, while the latter coincides with the rotation period derived 
from its light curve presumably caused by the presence of a spot 
or spots. Thus, this object is no longer a candidate showing pe- 
riodic TTVs (see Sec. 14.3b . We decided the leave it out from 
Table. [U 

KOI-895 A single-periodic companion around a solar-mass star 
with 1 1 .4 d TTV period. The stellar rotation, however creates 
a periodicity of 5.6 d, which is almost exactly half of the TTV 
period. This renders the main TTV period questionable. 
KOI-897 This planetary candidate is slighty larger than Jupiter 
and orbits a close-to solar mass host star with a period of 2.68 
days. There is a single period in its TTV Fourier-spectrum with 
82.0 days. 
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Table 1. Periodic TTV detections above 4-(T level. Periods in italic are ambiguous based on the analysis of the optical light curves of the candidates. 
The table contains 19 systems, 16 of which shows one or multiple significant TTV periodicities. 
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Fig. 5. A small part of the Kepler light curve of KOI-1152. The 
transits and eclipses are clearly seen. The eclipses occur at phase 
0.34 indicating an eccentric orbit. 



KOI-977 The host star of this object has a mass of 0.2 IM© and 
a radius of 16.48/?© with an effective temperature of 4204 K ac- 
cording to the Kepler Input Catalog. The nature of this star was 
corroborated by Muirhead et al. 2012 who clearly identified it 
as a giant based on near-infrared spectroscopic observations. The 
radius of the companion is correspondingly large, 63.45/?£, the 
largest in our sample. The mass of the host star is obviously too 
low, most probably a result of the procedure that was optimized 
to select main-sequence targets for the Kepler Mission. 
KOI-1003 This is monoperiodic TTV candidate with a long but 
highly significant period (277 days). 

KOI-1152 A larger-than- Jupiter companion orbits an early Ml 
dwarf (Muirhead et al. l20121 i. The star shows a 3 day rotation pe- 
riod. Besides the deep transits, the secondary minima (eclipses) 
are also easily seen. In addition, the orbit of the companion is 
obviously eccentric, since the eclipses are detected around 0.34 
phase instead of 0.50 (see Fig.|5]). The lagged secondary eclipse 
refers to a minimum excentricity of 0.26 (Eq. (2) in Dong et 
al l2012l) . The circularization time scale for this system is in the 
order of Tdrc ~ 120-270 Myr, assuming a planet mass of 1 Mj, 
a tidal quality factor of 10^ and other necessary system param- 
eters as derived in KIC and in Batalha et al. 2012 The shorter 
time scale assumes an excentricity of 0.26, the longer time scale 
is the e ~ Q limit. Because the circularization time scale is much 
less than the expected lifetime of an M, - 0.58Mq star, it is rea- 
sonable that the system is older than t„,c, and the excentricity 
is maintained by a dynamical process. Therefore, this system is 
a good candidate to have an outer companion that perturbs the 
orbit. 

KOI-1285 The optical light curve of this multi-frequency TTV 
candidate shows peaks close to the 53 and 70 d TTV periodici- 
ties, the other 3 long-period variations are not affected. 
KOI-1382 This system shows two very clear periodicites (with 
34.4 and 17.4 days), close to but not exactly at a 1 :2 period ratio. 
The rotation period of the I.IM© host star is 4.8 days. 
KOI-1452 The size of the companion of this star is twice that 
of the Jupiter. The system exhibits multiple TTV periodicities. 
The 74.4 and the 58.3 d peaks are close to peaks detected in 
the Fourier spectrum of the photometric light curve, therefore 
are ambiguous. The stellar rotation period is estimated to be 1.5 
days based on the Kepler light curve. The system survived all 
our criteria, and remains a strong multiperiodic TTV candidate. 



KOI-1540 If indeed a planet, this is one of the largest candidate 
in our vetted sample with almost three Jupiter radius, orbiting a 
solar-mass star One of the six significant periods (33 d) is found 
to be spurious from the light curve analysis. 
KOI-1543 KOI-1543 is a monoperiodic TTV candidate with a 
3.96 d orbital and a 97.00 d TTV period around a star of 1 .06M0 
mass and 0.86/?© radius. 

KOI-1546 This object features the richest spectrum of transit 
timing variations in our sample. Three frequencies (62.6, 357.9, 
74.3 d) are uncertain, because these appear in the light curve as 
well, another ten peaks reach the 4-cr level. 

4. Interpretation 

The observation of periodic TTV variations in hot Jupiters 
has been completely unexpected in the light of the loneliness 
paradigm, and would imply that one of the following three sce- 
narios acts in systems with hot Jupiters, and affects the transit 
times: 

- systems with transiting hot Jupiters can have non-transiting, 
massive close-in planets on highly inclined orbits; 

- hot Jupiters with significant TTVs host large moons; or 

- some other, still unidentified source of periodic TTV is in 
action. 

Here we discuss these possibilities in more details. 

4.1. Are there more planets in there? 

Hot Jupiters are not necessarily single. If that is the case, peri- 
odic TTVs may be tracers of perturbations from additional plan- 
ets in the system, which are not observed in transits because of 
the viewing geometry. A known case for a non-single hot Jupiter 
is in HAT-P-13, hosting a hot Jupiter and a distant massive planet 
on eccentric orbit, which is not observed in transits (Szabo et al. 
|20T0] |. In the Batalha et al. l20T2l catalog of 2321 planet can- 
didates, there are several examples for hot Jupiters in multiple 
systems, e.g. KOI-338, a system with two planets, both are hot 
Jupiters. KOI-94 consists of four transiting candidates, one hot 
Jupiter, two hot Neptunes and a hot super-Earth, KOI- 1241 hosts 
a hot Jupiter and a hot Neptune; etc. 

However, the distribution of hot Jupiter candidates with TTV 
is very different from the distribution of planets in multi sys- 
tems. In the left panel of Fig. |6] we plot the apparently sin- 
gle hot Jupiters with crosses and known multiple systems with 
stars (the left upper part of the region they occupy). We found 
that very few multiple systems exhibit planets above the line 
R/Re > 5 yjP/day. On the other hand, all but one TTV hot 
Jupiters reside above this line, while they follow the same distri- 
bution as non-TTV, (apparently) single hot Jupiters. One could 
argue that massive, non-transiting companions on an inclined or- 
bit can perturb the transiting planet, but in this case, the planets 
must be also close-in types to support the median detected TTV 
period of 50 days. These systems would look very differently 
from e.g. HAT-P-13. 

4.2. Are there exomoons in there? 

Another - however, still unconfirmed - explanation of TTVs is 
the presence of moons, pumping TTVs (Sartoretti and Schenider 
fT9991 Szabo et al. l2006l Simon et al. l2U07l Kipping [2009b . The 
expected rate of a TTV can be estimated if we assume that the 
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satellite has no sign in the light curves, only affects the tran- 
sit time of the planet. In this case, the TTV full amplitude is 
a sM sP{mM pY"^ , where a, and a is the semi-major axis of the 
satellite and the planet, M, and Mp is their mass, respectively, 
and P is the orbital period (Sartoretti and Schneider, 1999). For 
an order-of-magnitude estimate of the expected effects, we as- 
signed 1 My mass to the planet, 1 Me mass for the exomoon, 
and assumed a, - an, i.e. the moon orbits at the Hill-radius. 
With these assumptions, TTV caused by a moon can be calcu- 
lated as a function of the orbital period of the planet. 

In the right panel of Fig. |6] we plot the measured TTV am- 
plitudes for the hot Jupiters. Systems with multiple TTV peri- 
ods are plotted with different symbols, and they follow the same 
trend as the monoperiodic systems. The linear model overplots 
the expected TTV amplitude of a Jupiter + Earth system. We 
conclude that the distribution of TTV points above an orbital 
period of 2 days follows a similar trend as the model, but for 
shorter orbital periods, TTVs significantly exceed this predic- 
tion. However, TTV can be underestimated by a factor of 2-5 
based on the Sartoretti and Schneider (1999) assumptions espe- 
cially when the moons are quite large (see Szabo et al. 120061 
Simon et al. 120071) . and the detected TTVs may be explained 
with less massive exomoons. 

An argument against the exomoon scenario is that a large 
(> 10%) fraction of hot Jupiters should have distant and mas- 
sive moons to explain the general appearance of the distribution 
(right panel of Fig.|6l). Besides challenging the origin and preva- 
lence of such systems, the interpretation also faces the question 
of stability. Hot Jupiters with a moon suffer tidal forces from 
both the star and from the close-in, massive exomoon. As a con- 
sequence, the rotation rate of the planet evolves, and the orbital 
distance of the exomoon can swing very significantly, either to- 
ward spiral-in, or escape, or oscillation, or chaotic behaviour 
(Barnes & O'Brien 2002). Taken this argument into account, one 
could raise a reasonable doubt about the assumption of massive 
exomoons around a significant fraction of hot Jupiters. 

4.3. Other scenarios 

It is possible that no more objects are necessary to invoke inter- 
preting the observed periodic or quasi-periodic TTVs. Another 
possible explanation is that the activity of the host stars, i.e. the 
presence of stellar spots perturb the transit time measurements. 
An excellent example is Kepler-17b (former KOI-203) which 



we originally included in our sample. Desert et al.|20TT| demon- 
strated that both the light curve and the O-C of the mid-transit 
times of this object are modulated by the stellar rotation period 
and its first harmonic. We also detected the second harmonic in 
the Fourier spectrum, and therefore discarded this object from 
further investigations. After scrutinizing the Kepler light curve 
of our target stars, one single periodic target (KOI-883) and 
about a dozen frequencies in multi-periodic TTV Jupiters were 
also eliminated because of the suspected indirect effect of the 
stellar rotation and activity. We note in passing that young, active 
stars show rotational period less than 10 days, and more aged, 
quiet stars can exhibit rotation periods up to 30 days. Modulation 
periods longer than that are unlikely to be caused by activity and 
rotation. 

Conclusions 

In this work we selected planetary systems from the Kepler 
sample having only one known hot Jupiter candidates without 
any other planets. Hot Jupiters are thought to be solitary planets 
orbiting their host stars. We searched for periodic TTV varia- 
tions in the sample, a possible sign of exomoons. After elimi- 
nating false positives using the results of Santerne et al. 120121 
and systems that might show spurious periodic signals due ef- 
fect of inadequate sampling and contributions from stellar rota- 
tion, we ended up with 16 systems. Among these objects 9 show 
only one significant period (KOI-131, 186, 256, 882, 897, 1003, 
1 152, 1448, 1543), others show multiple periods (KOI-823, 977, 
1285, 1382, 1452, 1540, 1546). The number of significant pe- 
riods ranges from two (KOI-977, 1382) to ten (KOI-1546). We 
caution that the detected periods are lower limits, since the sam- 
pling is based on the planetary transits in the system, and this can 
undersample any short-period variations in the Fourier-space. 
With a median mass of 1 .03 M© our sample consisting of mainly 
main-sequence dwarf stars closely follows the distribution of 
host stellar masses of the Batalha et al. 120 121 list (median: 1.05 
Mq) with two low-mass stars (KOI-256, 1152) and the evolved 
KOI-977. 

One of our best candidates is KOI-1152, which is an Ml 
dwarf star hosting a giant planet on a 4.7 d orbit. The eclipses 
can be detected, and the orbit is obviously eccentric. Since the 
circularization time is much less than the expected age of the the 
host star, we hypothesize that the eccentricity and the observed 
TTV period is maintained by a third body in the system. 
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Concerning the false positives in our sample, Santerne et al. 
120121 recently estimated that the rate of false positives among 
short-period Jupiters in the Kepler sample can be as high as 
34.8% + 6.5%. In the frame of their radial velocity follow-up 
survey none of our candidates were confirmed as planets, but 
a few were rejected as we discussed in Sec. 13.11 Previously, 
Morton & Johnson|20TT]estimated a rate which is closer to 5-10 
%. Depending on which value is taken, the non-planetary can- 
didates may be between 1 and 6 in our sample, even if brown 
dwarfs are permitted. More follow-up observations are clearly 
needed to settle this problem. It is worth noting that 4 of our 
candidates qualify for a vetting flag of 2, seven of them 3 and 
the rest (with the highest KOI numbers) have insufficient follow- 
up observations, therefore have a vetting flag 4 in Borucki et al. 

i2uni 

We discussed three possible scenarios to explain the ob- 
served periodic TTVs. If additional planets are present in the 
systems of these seemingly solitary hot Jupiters, they present 
a new population of non-transiting, massive, close-in planets, 
that has not been seen earlier We emphasize that the TTV hot 
Jupiters follow the distribution of solitary hot Jupiters on the 
period-radius diagram. 

Although the existence of large (M>1 Me) exomoons is an 
attractive possibility, such a scenario needs to address the for- 
mation and stability problems of large exomoons around many 
hot Jupiters. We also eliminated the possibility of the influence 
of sampling and stellar activity on the observed transit timings. 

Longer observations and ground-based follow-ups can help 
in confirming the nature of the TTV-positive hot Jupiters. With 
the help of further Kepler short cadence observations, or even 
with accurate transit photometer with large telescopes applying 
much denser sampling than Kepler, sampling effects in long ca- 
dence data can be eliminated. With this technique we will be 
able to confirm that the observed TTVs are not due to effects 
of regularly spaced sampling with long exposures. Even Kepler 
LC data can help us confirm that the observed TTV is not due to 
stellar rotation and activity. That is because the appearance of the 
stellar surface evolves with the activity cycle, thus the effect will 
be different in different phases of the activity cycle. If the phys- 
ical nature of the observed hot Jupiter TTVs can be confirmed 
with appropriate significance, only dynamical explanations sur- 
vive: additional planets on exotic orbits, or the exomoons. 
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